Local Behavior of the First-Order Gradient Correction to the Thomas-Fermi Kinetic 

Energy Functional 
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The first order gradient correction to the Thomas-Fermi functional, proposed by Haq, Chattaraj 
and Deb (Chem. Phys. Lett. 81, 8031, 1984) has been studied by evaluating both the total kinetic 
energy and the local kinetic energy density. For testing the kinetic energy density we evaluate its 
deviation from an exact result through a quality factor, a parameter that reflects the quality of 
the functionals in a better way than their relative errors. The study is performed on two different 
systems: light atoms (up to Z — 18) and a noninteracting model of fermions confined in a Coulombic- 
type potential. It is found than this approximation gives very low relative errors and a better 
local behavior than any of the usual generalized gradient approximation semilocal kinetic density 
functionals. 

PACS numbers: 31.15.Ew, 31.15.Bs, 31.15.-p, 71.10.CA, 71.15.Mb 



I. INTRODUCTION 

Density Functional Theory, originally developed by 
Hohenberg and Kohn (HK), 1 replaces the main role of the 
wave function in the description of many electron systems 
and proves that the particle density n (r) determines all 
the properties of the system through the functional of the 
total energy E [n] . When E [n] is minimized, the density 
and the energy of the ground state can be found. The 
functional for the total energy of the interacting system is 
unknown, but Kohn and Sham (KS) 2 devised a method 
for the minimization of E [n] that uses the concept of non- 
interacting KS orbitals <pi (r) that yield the same density 
than the interacting system. The KS scheme allows to 
divide E [n] in four terms, and one of them is the kinetic 
energy density functional (KEDF) of the KS noninter- 
acting system, Tg [n] , that can be exactly evaluated in 
terms of the KS orbitals. 

In order to provide insight into the physics of many- 
body systems, we are interested in elucidating the for- 
mal properties of the KEDF. As explained below, the 
Thomas-Fermi (TF) functional is a local density approx- 
imation that has been corrected with other terms that 
usually depend on the variations of the density, normally 
including second-order gradient terms (see Ref. [1 for a 
general review on the results obtained with a number 
of semilocal kinetic functionals). In this paper we will 
study the kinetic energy functional that includes the first- 
order gradient correction to the TF functional proposed 
by Haq, Chattaraj and Deb4^>^ To assess the qual- 
ity of this first-order gradient correction we will evaluate 
both the total kinetic energy and the local behavior of the 
approximate kinetic energy density (KED). The study of 
the KED (a function that depends on the spatial coordi- 
nates) is made through a quality factor a, a parameter 
that reflects the quality of the functionals in a better way 



than their relative errors. 3 ' 9 ' 10 ! 11 ' 12 It may be expected 
that the functional which systematically gives the better 
KED would probably also yield another properties more 
accurately. 

In order to make the study of how this first-order gra- 
dient functional works for localized systems we test two 
simple systems: a noninteracting model (fermions con- 
fined in a Coulombic-type potential) and an interacting 
electron system (light atoms, from He to Ar). 



A. Kinetic Energy Functionals 

As commented, the KS scheme allows to define the 
KEDF of the KS noninteracting system, T$ [n] , for any 
many particle system. But, following the original HK 
idea, many attempts to write the kinetic energy as a 
functional, using exclusively the electron density, have 
been made. The very first attempt was a statistical 
method - the Thomas-Fermi (TF) functional, proposed 
by Thomas^ 3 - and Fermiii - that has the form: 

Ttf M = C TF J drn(r) 5 / 3 , C TF = ^(3^ 2 ) 2 / 3 . 

(1) 

TF is exact for a free electron gas or an uniform system, 
but this local density approximation yields usually big 
relative errors -around 10%- for systems such as atoms 
and molecules. 

Another well-known functional, due to von Weizsacker 
(vW)^, is exact for any system described by a single 
orbital, 



T vW [n] 



1 



dr 



|Vn(r)p 
n(r) 



(2) 



The TF functional is usually corrected by function- 
als that depend on the variations of the density, as in 



2 



the second-order Gradient Expansion Approximation) 16 
(GEA2, where the vW term is weighted by 1/9) 



T GE A2 [n] = Ctf Jdrn{r) 5 / 3 + 



72 / dT n(r) ' [6) 



that usually yields small relative errors -usually lower 
than 1%— when applied to localized systems using good 
densities (i. e., those obtained with accurate methods 
as the Hartree-Fock or the KS ones), although it gives 
a wrong local behavior, misplacing the KED about 
18%4 Another functional that has obtained good re- 
sults for these type of systems (when atomic systems 
and homonuclear molecules are variationally solved) 17 > 18 
is the TF±vW one, 



T, 



TF±vW l n \ 



-T v w [n] 



T T f [n 
Ctf J dr n(r) 



40 J n(r) 



1(4) 



But if we use these good electron densities, TF-^vW usu- 
ally yields bigger relative errors than GEA2 and usually 
overestimates the total kinetic energy by more than 6%. 

Beyond the GEA2 and TF±vW, a number of other 
semilocal functionals within the Generalized Gradient 
Approximation (GGA) can be proposed using the gen- 
eral form 



Tgga [n] = Ctf J dr n(rf/ 3 F [ 



(5) 



where F [s (r)] is an enhancement factor that corrects the 
TF functional via the adimensional reduced gradient, 



S (r) 



|V"(r)| 
n(r)V3 ■ 



that takes into account the inhomogeneities of the sys- 
tem. The GGA functionals include both the TF func- 
tional (with F [s (r)] = 1) and the GEA2 and the TF±vW 
ones, as well as many other semilocal functionals. The 
GGA functionals generally depend on s(r) 2 , and their 
mathematical form and their results are cleary related 
to the GEA2; only some of the GGAs functionals yield 
smaller relative errors than the GEA2 but none of them 
give better KED than the TF functional 4 

It is also important to know that, due to symmetry 
considerations in general systems, these functionals de- 
pend only on the absolute value of the gradient of the 
density (or its squared value) and never take into account 
the directionality of the density variations. 



B. First-Order Gradient Correction 

For localized systems a different functional has been 
proposed by Haq, Chattaraj and Deb (HDC) -. This func- 
tional is an interesting exception, because it depends on 



the gradient of the electron density and not on its squared 
gradient, and also depends on the position vector r (dif- 
ferently to the other GGAs). 

The form of the proposed functional is: 



T HCD -i[n{v)] = Ctf J dr n(r) 5 / 3 - ^ J dr 



T TF [n (r)]+T r _![n (r) 



rVn(r) 
(6) 



i. e., the sum of the TF functional and an r-dependent 
correction T r _i [n (r)]. Note that the 1/40 weight resem- 
bles the same prefactor that appeared in the TF-ivW 
functional and that the gradient of the electron density 
is included in the r-dependent term, so this correction is 
usually labeled as first- order gradient correction (instead 
of the second-order corrections previously presented) . 
Making use of the relation 

, rVn(r) f , n(r) 
dr — = — / dr — — , 



a different first-order term is obtained, and a correspond- 
ing functional can be proposed, 

The functional T HCD -2 [n (r)] = T TF [n (r)]+T r _ 2 [n (r)] 
yields the same total kinetic energy but has a different 
local behavior. 

It is also known that T r _i[n(r)] without the 1/40 
prefactor is a lower bouncU^ to the Fisher information 
function^^I F \n{r)], 

, rVn(r) f , |Vn(r)| 2 r r . 

dr ^ < / dr 1 V ;l = I F [n r . (7) 

r z J nyr) 

On the other hand, the Fisher information function is 
directly related related to the vW functional, 



T vW [n] = h F [n (r)] 



(8) 



so |T„vkM i s an upper bound to T r _i[n.(r)] and 
T TF ± vW [n (r)] becomes an upper bound to the HCD 
functional, 

Thcd {n(r)}<T TF i vW [n (r)] . 

For all these results, the functional Thcd [n (r)] could be 
a promising functional yielding small relative errors. 

Note that T r _2 [n (r)] is the average of r~ 2 over the 
system, (r~ 2 ) = J dr^^-, with a 1/40 prefactor. This 
quantity is related to exact bounds for the von Weizsacker 
functional, being particularly interesting the following 
one: 



l(r- 2 )<IW [n(r)]<l(r- 2 ) 



(9) 



The lower bound is of variational origin while the upper 
bound is based on the complete monotonicity hypothesis. 
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(see a general discussion on this topic in Ref. (12 
references therein). 



and 



We must note that the choice of 1/40 as a prefactor has 
no theoretical justification and was chosen by trial and 
error. But it is possible to make a qualitative exercise to 
make a guess of the approximate value of this prefactor. 
We can argue using the previous inequality that the ki- 
netic energy obtained with the vW functional should be 
close to the mean of the two bounds, i.e. 



T vW [n (r)] . 



On the other hand, it is well known that the TF func- 
tional can be quite accurately corrected by adding a 
weighted vW term (e.g. the GEA2 functional usually 
yields errors lower than 1%). And by replacing the vW 
functional by A (r ~ 2 ), we add to the TF functional a 
correction that is i ( r ~ 2 )- This value 1/48 is quite close 
to the proposed 1/40 (a difference of 17%, for a correction 
that represents about 8% of the total kinetic energy). 



of valid KEDs using a linear combination of tg and tjJ , 



V ( r ) 



= (l-a)4(r)+ad/(r) 
1 N 

= (l-a)2£|V^(r)| 2 

i=i 

N ( 1 

(-oPi ( r ) y2 ^ W 



(11) 



(12) 



( 10 ) (all of them are valid KED, see Ref. (0)). So, t L s ' a {r) 
has a first term tg (r) given by the orbitals plus a term 
proportional to the laplacian of the electron density (the 
laplacian integrates to zero in the whole space). The real 
parameter a can be arbitrarily varied and the family of 
KEDs described by tg' a (r) includes tg (r) (when a = 0) 
and tjJ (r) (for a = 1). 

When an approximate KEFD is used, each value of a 
defines a quality factor 



cr(a) 



fdr t^(r)-t^ nc (r) 



T s [n] 



(13) 



The Relative Errors and the Local Kinetic 
Behavior 



As pointed out before, the present paper is devoted to 
the study of the two functionals proposed in Sect. pB[) . 
Thcd [n (r)L with an analysis of both their total kinetic 
energy and their KED. We study two different systems 
with spherical symmetric densities. Firstly, we use the 
electron densities of light atoms from Z = 2 to Z = 18. 
As a second system, we will study the noninteracting 
fermion system confined in a Coulomb- type potential pre- 
sented in Ref. (12). In both cases, we have worked with 
the analytical expressions of the single-particle orbitals in 
order to evaluate not only the relative errors obtained by 
the approximate KEDFs, but also compare their KED 
to the exact KS results (that can be directly obtained 
from these orbitals). To make that comparison we have 
introduced ^ 10 ' 11 ! 12 a quantity a that takes into account 
the local differences, focuses on the error cancellations, 
and reflects the fact that the KED is not univocally de- 
fined. 

Let's do a brief review of the method. For a given set of 
N orbitals (pi(r) we can use the usual orbital-based KED, 
tg (r) = i \ V<Pi( r )\ 2 ' or another common definition 



of the KED. 



3Ef=i^*(r)V 2 ^(r). But for 



this set of N orbitals we can construct an infinite number 



being t s unc (r) the approximate KED of the kinetic func- 
tional, a represents a summation of the absolute value 
of the differences between KEDs, so it measures how the 
approximated KED is misplaced with respect to the dis- 
tribution tg' a (r). Note that a is zero if the approximate 
and exact distributions of KED are identical, and when 
the larger the value of a is the bigger the differences be- 
tween the two are. For each approximate functional we 
will obtain a set of values of a (a) and in order to com- 
pare the KEDs we choose the best value of a for each 
functional, i. e. the value of a m i n that yields the min- 
imum value of a (a); the KED has its closest ig(r) by 
evaluating Eq. (fT2"|) for this value a m i n . For details on 
this measure a, the behavior of GGA functionals and the 
non uniqueness in the definition of the exact KED, see 
Ref. O) and references therein. 



II. RESULTS AND DISCUSSION 

In order to assess the quality of the two HCD, r- 
dependent functionals, introduced in Sec. II Bl (HCD-1 
and HCD-2), we are going to compare both the relative 
errors in the kinetic energies and the quality factor a of 
the KED obtained with these functionals to the values of 
the same magnitudes obtained with the TF, GEA2 and 
TF±vW functionals. 



We have applied the functionals to two different sys- 
tems. Firstly we study the atoms from He to Ar (from 
Z = 2 to Z = 18). We will approximate the atomic 
orbitals by orthogonalized Slater orbitals, because they 
do not add any spurious oscillations in the laplacian of 



a 
o 

a. 



0.10 



0.05 



-0.05 



-0.10 



-0.15 




0.24 -i 



0.22 



0.20' 



0.18- 



0.14 



0.12' 




~~ i — 1 — i — ' — i — ' — i — 1 — r 
10 12 14 16 18 



FIG. 1: (Color online) Relative errors of the TF, GEA2, 
TFjrvW and the first order correction HCD functionals for 
atoms from He to Ar. 



FIG. 2: (Color online) Quality factors a obtained with the 
TF, GEA2, TF±vW and the first order correction HCD func- 
tionals for atoms from He to Ar. 



the density and both the correct cusp behavior and the 
density decay for r — > oo can be achieved (the Gaussian- 
type basis sets have well known deficiencies, they cannot 
reproduce the cusps of electronic density and electronic 
orbitals on nuclei and cannot provide the correct asymp- 
totic behavior at long distances from the nuclei). Being 
the orbitals known, the local behavior of the KED can 
be exactly studied. The values of the exponents for the 
Slater orbitals 2 — >2£ and the corresponding exact total ki- 
netic energy for the atoms used in this paper are shown 
in Table [J 

Fig-Q]shows the values of the relative errors of the total 
kinetic energy yielded by the functionals for the atoms, 
as a function of the atomic number Z. The best results 
are obtained for the GEA2 and the two HCD functionals, 
with quite close errors. The TF functional gives a nega- 
tive relative error, whereas the TF^vW shows a similar 
trend, but with positive errors and smaller absolute val- 
ues. In all cases, the relative errors for all theses KEDFs 
are almost constant for Z > 10. 

For the atoms, the value of a is smaller the bigger Z 
is, as shown in Fig. [21 where the values of the quality 
factor are plotted as a function of Z. The results clearly 
show that the HCD-2 functional give the smallest values 
of a, reflecting the better local behavior of this semilocal 
functional. We have found that the GGA corrections to 
the TF functional improve the TF results for the total 
kinetic energies, but they always give worse local KEDs, 
an unexpected failure that shows that GGA functionals 
are unable to improve the local pathologies of the TF 
functional.- This failure is overcome by the HCD-2 form 
of the first-order correction, whereas the HCD-1 form 



improves the total kinetic energies, but not the value of 
a, as can be seen in Figs. [T]and[2j 

The second test system is constructed with a set of 
nonintcracting fcrmions that are confined in a Coulom- 
bic potential^ 2 - In this case, we can also use analyt- 
ical single-particle orbitals, those of a Coulombic-type 
potential. We will use spherically symmetric densities 
for this iV-fermion system built up with mixed states 
and we will describe the degenerate iV-fermion states 
that arise from the distinct values of the occupancy of 
the I single-fermion states. These degenerate 7V-fermion 
states have a different orbital angular momentum for 
each configuration, with a corresponding different den- 
sity n (r) As the HK theorems hold not only for 
densities that correspond to a ground-state density for 
some external potential^ but also for a weighted average 
of ground-state electron densities corresponding to that 
potential ) 27 ' 28 ^ 9 ' 30 we can directly discuss the densities 
corresponding to a degenerate ground-statei 30 ' 31 i 32 ' 33 

For the sake of simplicity we will only present the de- 
tailed results for the noninteracting systems with a num- 
ber of fermions N between 1 and 10. In TableHTlthe occu- 
pation numbers of each of the subshells of single-fermion 
orbitals Is, 2s and 2p, and the value of 53<=i \Jh(U + 1) 
for each of the A^-fermions configuration are shown. As 
commented in Ref. (12), y/hiji + 1) can be intu- 

itively related to an angular momentum of the configu- 
ration, although it is not the expected value of the angu- 
lar momentum operator but the sum of the eigenvalues 
of \[T? for each independent single-fermion orbital. In 
the table we also show the relative errors e obtained with 
the KEDFs and the values of a m i n that give the closest 
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TABLE I: Values of the Slater orbital exponents and the corresponding exact total kinetic energy for the atoms studied in the 
paper. 



Atom 


He 


Li 


Be 


B 


C 


N 


O 


F 


Ne 


Is 


1.6875 


2.6906 


3.6848 


4.6795 


5.6727 


6.6651 


7.6579 


8.6501 


9.6421 


2s 




0.6396 


0.9560 


1.2881 


1.6083 


1.9237 


2.2458 


2.5638 


2.8792 


2p 








1.2107 


1.5679 


1.9170 


2.2266 


2.5500 


2.8792 


Exact KE 


2.8477 


7.4421 


14.556 


24.498 


37.620 


54.281 


74.538 


98.996 


127.80 


Atom 




Na 


Mg 


Al 


Si 


P 


S 


CI 


Ar 


Is 




10.6259 


11.6089 


12.591 


13.5745 


14.5578 


15.5409 


16.5239 


17.5075 


2s 




3.2857 


3.6960 


4.1068 


4.5100 


4.9125 


5.3144 


5.7152 


6.1152 


2p 




3.4009 


3.9129 


4.4817 


4.9725 


5.4806 


5.9885 


6.4966 


7.0041 


3s 




0.8358 


1.1025 


1.3724 


1.6344 


1.8806 


2.1223 


2.3561 


2.5856 


3p 








1.3552 


1.4284 


1.6288 


1.8273 


2.0387 


2.2547 


Exact KE 




161.14 


198.82 


242.86 


289,19 


340.88 


397.34 


458.88 


525.59 



TABLE II: Relative errors e and optimized values of the quality factor a m in for all the configurations of the noninteracting 
iV-fermion system from N = 1 to 10. The occupation numbers p of each orbital is shown, as well as the values of Yli=i h(h + l)- 





Orbital occupation 


TF 


TF^vW 


GEA2 


HCD-1 


HCD-2 


N 


Pis P2s P2p J2i=l + 1) 


min C 


O'min £ 


& min c 


O'rnin £ 


min c 


1 


10 


0.422 -0.422 


0.222 -0.222 


0.311 -0.311 


0.322 -0.322 


0.330 -0.322 


2 


2 


0.166 -0.082 


0.225 0.118 


0.187 0.029 


0.171 0.018 


0.137 0.018 


3 
3 


2 10 
2 1 2 


0.179 -0.107 
0.198 -0.126 


0.236 0.073 
0.244 0.041 


0.200 -0.007 
0.215 -0.033 


0.182 -0.013 
0.201 -0.036 


0.148 -0.013 
0.170 -0.036 


4 
4 
4 


2 2 
2 11 2 
2 2 4 


0.180 -0.114 
0.196 -0.132 
0.209 -0.145 


0.233 0.057 
0.238 0.026 
0.241 0.003 


0.199 -0.019 
0.211 -0.044 
0.220 -0.062 


0.182 -0.024 
0.197 -0.046 
0.210 -0.061 


0.149 -0.024 
0.166 -0.046 
0.182 -0.061 


5 
5 
5 


2 2 1 2 
2 12 4 
2 3 6 


0.187 -0.127 
0.197 -0.139 
0.206 -0.147 


0.227 0.026 
0.228 0.003 
0.227 -0.013 


0.201 -0.042 
0.207 -0.060 
0.212 -0.073 


0.187 -0.043 
0.197 -0.059 
0.205 -0.070 


0.158 -0.043 
0.169 -0.059 
0.180 -0.070 


6 
6 
6 


2 2 2 4 
2 13 6 
2 4 8 


0.182 -0.127 
0.188 -0.135 
0.193 -0.139 


0.213 0.012 
0.209 -0.004 
0.206 -0.016 


0.192 -0.049 
0.194 -0.062 
0.195 -0.070 


0.181 -0.049 
0.187 -0.060 
0.191 -0.067 


0.155 -0.049 
0.161 -0.060 
0.168 -0.067 


7 
7 
7 


2 2 3 6 
2 14 8 
2 5 10 


0.169 -0.117 
0.172 -0.122 
0.173 -0.123 


0.193 0.011 
0.186 -0.001 
0.180 -0.008 


0.176 -0.046 
0.174 -0.055 
0.173 -0.059 


0.168 -0.044 
0.170 -0.052 
0.171 -0.055 


0.143 -0.044 
0.147 -0.052 
0.150 -0.055 


8 
8 
8 


2 2 4 8 
2 15 10 
2 6 12 


0.151 -0.101 
0.150 -0.103 
0.149 -0.102 


0.178 0.019 
0.158 0.010 
0.150 0.006 


0.156 -0.035 
0.150 -0.040 
0.146 -0.042 


0.148 -0.032 
0.147 -0.037 
0.146 -0.038 


0.126 -0.032 
0.126 -0.037 
0.127 -0.038 


9 
9 


2 2 5 10 
2 16 12 


0.131 -0.080 
0.125 -0.080 


0.169 0.032 
0.143 0.027 


0.145 -0.018 
0.127 -0.020 


0.132 -0.015 
0.122 -0.017 


0.110 -0.015 
0.102 -0.017 


10 


2 2 6 12 


0.123 -0.056 


0.166 -0.051 


0.141 0.004 


0.126 0.007 


0.106 0.007 



tg' a (r) in Eq. (fT2|) for each functional. For the TF func- 
tional, the relative errors get smaller when the number 
of fermions increases. Showing much smaller absolute er- 
rors (probably due to a compensation of the TF errors by 
the opposite trend of the vW functional: its errors grow 
with the number of fermions the TF±vW functional 
yields a similar behavior, except that the errors slightly 
increase for N > 8. The GEA2 and both HCD func- 
tional present no clear trends, because they appear to 
work better for those iV-fermion configurations that are 
close to the closed-shell. And for all the functionals but 
the TF-ivW functional the relative error is smaller for 



configurations with smaller values of $^ i=1 y/h{h + !)• 
Moreover, note that the errors are much smaller for the 
TF-|vW functional than for the GEA2 one, an opposite 



trend that we found for light atoms. 

As explained in Section [I CI u m in is the amount of ap- 
proximate KED misplaced with respect to the distribu- 
tion tg' a (r) after minimizing with respect to the param- 
eter a — see Eq. (fT5)l . The values of the quality factor 
of Table [TX| show that a m in decreases with the number of 
fermions for the TF functional. Even the relative errors 
of the kinetic energy are much smaller than for the TF 
case, the TF-jrvW functional yields values of a that are 
bigger than those obtained for the TF functional, follow- 
ing the general results shown by the GGA functionals. 3 
Again, the local compensation of the TF behavior and 
the opposite trend of the vW functional^ cause the re- 
sult for the TF-jrvW functional. These results for cr m , n , 
and the previous mentioned relative errors, arise from 
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FIG. 3: (Color online) Relative errors obtained with the TF, 
GEA2, TF|vW and the first order correction HCD function- 
al in terms of 52j=i y/U{U + 1), for all the 64 configurations 
of the confined noninteracting system with 34 fermions. Note 
that these results do not discriminate configurations with dif- 
ferent values of p4 S , as shown in the results for a of Fig.@ 



FIG. 4: (Color online) Quality factor a obtained with the 
TF, GEA2, TF±vW and the first order correction HCD func- 
tional in terms of Yli=i ^h(h + 1), for the 64 configurations 
of the noninteracting 34-fermion system. Note the different 
results the two HCD functionals yield. 



the known properties of the vW functional: it is a lower 
bound to both the total KE and the local KED. On the 
other hand, the general trend for the GEA2 and both 
first-order HCD functionals indicates that the quality fac- 
tor decreases with the number of fermions, with quite 
similar results for the three approximations. 

Focusing on the dependence on the angular momentum 
J2iLi \/h(k + 1) of the 7V-fermion configurations we see 
that the local description of the TF functional improves 
for smaller angular momentum (i. e. the quality fac- 
tors decrease with the angular momentum). The TF-^vW 
functional shows the opposite trend. On the other hand, 
the values of a for the GEA2 and the two HCD function- 
als decrease with the angular momentum when the shell 
is less than half-filled and increases when it is more than 
half-filled. 

In order to make some other comparisons on the de- 
pendence on the angular momentum, let us discuss the 
case with N = 34 fermions, that has 64 different config- 
urations, all degenerate with the same energy. As com- 
mented before, these degenerate A^-fermion systems have 
spherically symmetric densities, but due to the different 
occupancy of the I single-fermion states they give differ- 
ent results for the kinetic energies and for the values of 
a. Now we will focus on how the results obtained with 
each of the different KEDFs depend on J2iLi \/h{k + 1)- 
Regarding the relative errors, the HCD r-dependent func- 
tionals obviously give the same value since, once in- 
tegrated, both yield the same kinetic energy. Fig. [3] 



shows the values of the relative errors for the function- 
als for the noninteracting fermion system, as a function 

°f J2iLi \/h(k + I). In this case, the HCD r-dependent 
functional gives the best results, better than the GEA2 
ones, whereas the TF^vW functional improves the values 
obtained with the TF one, getting again positive values. 
Note that for this system all the values of the relative 
errors are much smaller than for the atomic system. 

Fig. 0] depicts a as a function of Yli=i \/h(h + !)• hi 
this case, the TF and HCD-1 functionals give very close 
results for a, whereas the GEA2 functional yields worse 
values than the TF and the results of the TF^vW func- 
tional clearly go bigger. To better understand the results, 
let us analyze how the fermions distribute among all the 
single-fermion orbitals. For A" = 34, the fermions oc- 
cupy the lowest-energy orbitals, with six fermions in the 
fourth shell. Now we put our attention on the configu- 
rations that have a fixed value of X^iLi \/h(h + 1), i- e. 
those results with the same value for the abscissas in the 
figure. Using the arguments presented in Ref. (fl2l). for 
a fixed value of \/h(h + 1) the configurations with 

ri4 S = 2 have more higli-^ single-fermion orbitals occu- 
pied than the configurations with m s — 1; and when 
m s = we will have our six fermions occupying orbitals 
with the lowest-Zi possible values. Fig. 2]then shows that 
the values of a for the distinct configurations are clearly 
ordered in terms of the different values of p^ s . This is 
not the case for the relative error of the energies, as can 
be seen in Fig. [31 where no differences are observed. This 
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effect was afso previousfy studied in Ref. (fl2l). 

III. CONCLUSIONS 

The first-order correction HCD functionals improve the 
semilocal functionals (GEA2 and TF-ivW) we have used 
as benchmarks in the test systems we have studied. The 
relative errors are quite similar to the GEA2 functional 
for atoms, whereas they are slightly smaller for the non- 
interacting confined fermion system. We get accurate 
total kinetic energies, but when the local KEDs behavior 
is measured through a, only the HCD-2 functional clearly 
gives better values than the TF functional, whereas the 
HCD-1 functional gets quite close results to those of the 
TF one. These results have revealed the differences in the 
local behavior of both HCD functionals. On the other 
hand, Figs. [2] and Q] show that the values of a are clearly 
smaller for the noninteracting fermion system than for 
the atomic system, reflecting the importance of the in- 
teraction on the KED. 

As was also found in the previous paper,— we have 
shown that for a specific number of fermions the results 
for the KEDs generally are different for each of the de- 
generate states of the noninteracting ./V-fermion system 
with spherical symmetric densities. Depicting their val- 



ues of a versus the values of 2i=i \/h(k + 1) of each 
specific configuration some trends appear: for a given 
Si=i \/k{k + 1), those configurations whose occupied 
orbitals have the lowest values of k compatible this value 
°fJ2iLi y/hQi + 1) yield the smaller values for a. As a 
consequence, the KEDFs discussed in the paper give a 
better description of the KED when the lower-ij single- 
fermion orbitals have the maximum occupancy. These 
KEDFs give a better description of the KED when the 
lower-Zi single-fermion orbitals have the maximum occu- 
pancy. 

An open question is to understand the reasons why one 
of these first-order correction HCD functionals gives bet- 
ter local behavior than any other KEDF in the literature. 
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